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We consider nucleon-nucleon (A'A') interactions from chiral effective field theory. In this work we restrict 
ourselves to the elastic A'A' scattering. We apply the N/D method to calculate the A'A' partial waves taking as 
input the one-pion exchange discontinuity along the left-hand cut. This discontinuity is amenable to a chiral 
power counting as discussed by Lacour, Oiler and MeiBner [Ann. Phys. (NY) 326, 241 (2011)], with one-pion 
exchange as its leading order contribution. The resulting linear integral equation for a partial wave with orbital 
angular momentum i" ^ 2 is solved in the presence of i" — 1 constraints, so as to guarantee the right behavior 
of the D and higher partial waves near threshold. The calculated A'A' partial waves are based on dispersion 
relations and are independent of regulator. This method can also be applied to higher orders in the calculation 
of the discontinuity along the left-hand cut and extended to triplet coupled partial waves. 
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I. INTRODUCTION 

The nucleon-nucleon (A^A^) interaction is a basic process 
whose understanding is necessary for the study of nuclear 
structure, nuclear reactions, nuclear matter, neutron stars, etc. 
[1-5]. Since the early 1990s [6-8] the low energy effective 
field theory (EFT) of QCD, chiral perturbation theory (ChPT), 
has been applied to AW scattering in a large number of studies 
[6-18]. A sophisticated stage has been reached where the A^A^ 
potential is calculated in ChPT up to N^'LO [16, 18]. 

However, as the A'A^ interaction is nonperturbative, the 
chiral A^A^ potential must be iterated. It was proposed 
by Weinberg in his seminal papers [6, 7] to solve a 
Lippmann-Schwinger equation. Since the chiral potential is 
singular at the origin a regularization method, typically a 
three-momentum cut-off A [9, 16, 18], should be introduced 
for solving the Lippmann-Schwinger equation. Despite the 
great phenomenological success achieved by the AW chiral 
potentials in describing A^A^ scattering data [16, 18], a remnant 
cut-off dependence is left. It was shown in the literature 
[19-23] that the chiral counterterms introduced in the AW 
potential following naive dimensional analyses are not enough 
to renormalize the resulting AW scattering amplitude. In 
Ref. [19] one counterterm is promoted from higher to lower 
orders in the partial waves ^Pq, ^P2 and ^D2 due to the 
attractive 1/r^ tensor force from one-pion exchange (OPE). 
As a consequence, stable results, independent of cut-off in 
the limit A — > oo, are achieved when the LO OPE potential 
is employed. Partial wave amplitudes with larger orbital 
angular momentum £, £ ^ 3, can be calculated in Born 
approximation with sufficient accuracy [12, 24]. Then they 
do not pose any problem for renormalization, making use 
of standard perturbative techniques. It is also argued in the 
same Ref. [19] that higher orders terms in the chiral AW 
potential could be treated perturbatively. Reference [19] was 
extended along these lines to subleading two-pion exchange 



(TPE) in Ref. [25]. The promotion of higher orders to lower 
ones due to nonperturbative renormalizability is studied in 
detail in Ref. [24] by making use of the regularization group 
equations. (See also Refs. [20, 21, 26] for a coordinate 
space renormalization by imposing appropriate boundary 
conditions.) On the other hand, Refs. [27, 28], following the 
philosophy of Refs. [29, 30], stress that the cut-off A should 
not be taken beyond the breakdown scale of the EFT, typically 
around 1 GeV. It is argued that if this is done the power 
counting associated with the chiral EFT is lost. 

We employ here the N/D method [31] for studying A^A^ 
interactions in the elastic case. A linear integral equation 
then results for determining the NN partial waves. The input 
is given by the discontinuity of the partial wave along the 
left-hand cut (LHC), which is due to multipion exchanges, 
the lightest one being OPE. The well-known behavior of 
a partial wave near threshold, that vanishes like |pp, 
with IpI the center of mass (cm.) three-momentum, is not 
automatically fulfilled in the N/D method for £ ^ 2 [32-34]. 
Then, the N/D method must be solved in the presence of 
£ — I constraints. These are satisfied by introducing £ — I 
Castillejo-Dalitz-Dyson (CDD) poles [35], as we discuss 
below. The extension of the work to the inelastic partial waves 
will be considered separately [36], due to the intrinsic specific 
difficulties that occur in that case. The latter do not appear 
for the elastic partial waves which then allows us to focus in 
those aspects that directly concern the application of the N/D 
method to A^A^ scattering. 

After this Introduction we discuss the N/D method for 
calculating the A^A' elastic partial waves in Sec. II. Special 
attention is paid to derive the constraints needed to meet the 
threshold behavior for a partial wave with £^2. The inclusion 
of CDD poles for satisfying these constraints is a novelty 
in the literature. The results are discussed in Sec. Ill and 
conclusions are given in Sec. IV. The Appendix discusses the 
numerical method employed to solve the integral equation. 
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II. APPLICATION OF THE N/D METHOD TO 
ELASTIC A'A' PARTIAL WAVES 

A. A'A' partial waves cuts 

We consider the A^A^ interaction process 

A^(pi;(Jiai)A^(p2;(T2o:2)^A^(Pi;c7i'aOMP2;(^2«2) 
whose scattering amplitude in the cm. frame is indicated by 

{^\a[a[a2aJi\TdMi,a\aia2a2) . 

Here the initial momentum is p = |p|z, taken along the z 
axis, and the final one is p'. Its decomposition in partial 
waves is discussed in Appendix A of Ref. [1], to which we 
refer for further details. We denote a AW partial wave by 
7//5(^'>^; IpP). being (! the final orbital angular momentum 
and i the initial one. The labels J, S, and / stand for the 
total angular momentum, spin and isospin of the reaction, 
respectively: 

Tjis{i',t,\p\^) 

X {m's'-^sj,\e'SJ) {a[ 0^/3 1 Ti T2/) («! a2h I Tl Tit) 

^dp'{p\a[a[a^a!,\Td\\p\z,aiaia2a2)Yf?'ip'r . (1) 

In this equation, the Clebsch-Gordan coefficients for the 
couplings of two angular momentum y'l, J2 to 73 are indicated 
by {inim2ino,\jiJ2JT,), with oti, 1112 and OT3 the corresponding 
third components. 

A A^A^ partial wave amplitude has two cuts [37], the 
right-hand cut (RHC) for < p^ < 00, due to unitarity, and the 
LHC for — cxD < p^ < L with L = — m^/4, because of crossed 
channel dynamics. Both cuts are depicted in Fig. 1. The upper 
limit for the latter is given by OPE, as the pion is the lightest 
particle that can be exchanged in the t and u channels. Because 
of unitarity a partial wave amplitude satisfies in the cm. 
frame, above the elastic threshold and below pion production, 
the relation 




FIG. 1. The thick lines correspond to the RHC and LHC, in order 
from top to bottom. In the same figure the integration contours Cj 
and Cjj for evaluating Dj(>s{A) and Nj[s{A), respectively, are shown. 
One has to take the limit e — > 0+ 



to that of the reducible parts and it becomes a perturbation. 
At LO, according to the counting developed in Refs. [1, 4] 
(that for two-nucleon irreducible diagrams coincides with the 
standard chiral counting [6-8]), ImTjjs along the LHC is due 
to OPE. 
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with m the mass of the nucleon. In our normalization, the 
S'-matrix is given by 



Sjis — 



.'"IpL 



As shown in Ref. [1] one can calculate perturbatively 
in ChPT ImTjjs along the LHC, since this imaginary part 
is due to multipion exchanges. The infrared enhancements 
associated with the RHC are absent in the discontinuity along 
the LHC because, according to Cutkosky's theorem [38, 39], 
it implies to put on-shell pionic lines. Within loops the pion 
poles are picked up, making that the energy along nucleon 
propagators now is of ff{p), instead of a nucleon kinetic 
energy. In this way, the order of the diagram rises compared 



B. S and P waves (^ = 0, 1) 

In the following we take the elastic case for which i' = 
i = J (except for the ^Pq.) The N/D method [31] rests on the 
separation between the RHC and LHC. In this way, a partial 
wave Tj(s{A) is written as' 



TjesiA) = 



DjesiA) 



(3) 



The function Nj(s{A) has only LHC while Djis{A) has only 
RHC. Taking into account elastic unitarity, Eq. (2), one can 



' We replace the subscript / by f when denoting a partial wave. The former 
can be deduced from i and S. 
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ImDjtsiA) = -NMsiA)''-^ , A>0. (4) 

An 

Along the LHC one also has from Eq. (3) 

lmNm{A) = Dm (A) Im7>«(A) , A < -ml/4 . (5) 

We first write down a dispersion relation (DR) for Dji;s{A) 
and Nj(s{A) taking as integration contours C/ and C// in 
Fig. 1, respectively. The integration along the circle at infinity 
vanishes, if necessary, by taking the sufficient number of 
subtractions. For the case of a once-subtracted DR the 
following expressions result 



DjtsiA) = 1 - 



A-D 



NjtsiA)^Nm{D) + 



dq-" 



A-D 



p{q'-)Nj,s{q^) 
{q^-A){q^-D) ' 



Ak 



. Ms{k')Djts{k'-) 
{k^-A){k'--D) 



(6) 



(7) 



Here we have indicated by p{A) — m^/A/4-K for A > 
and Aj(>s{A) = lmTji>s{A) foi A < L — —ml/4. For physical 
scattering A ^ A + ie . The subtraction constant in Djfs{A) 
has been fixed to 1 because, according to Eq. (3), only the 
ratio between NjtsiA) and Dj(s{A) matters for determining 
Tjes{A)- In this way, one has the freedom to fix the value of 
Djes{A) at one point. 

Asymptotically, for p' -^ — oo, OPE tends to constant, 
so that, according to the Sugawara and Kanazawa theorem 
[32, 40] one subtraction is necessary for the DR of A'jfs(A), 
even though Aji:s{A) — > 1/A in the case of OPE. On general 
grounds, a partial wave amplitude is bound for A ^ +cxd by 
a constant because of unitarity and the same theorem then 
requires that at least one subtraction is necessary. 

An integral equation for Djes{A) results by inserting the DR 
for NjisiA), Eq. (7), into Eq. (6): 



Djes{A)^l-Nj(s{D) 



A-D r+ 



A-D 



dq' 



Piq') 



K Jo 



Ak 



dq^ 



Piq") 



{q2-A){q2-D) 



Ajes{k^)Djesik^ 



n^ Jo ^ q^-Aj-^ {k^- - q^ - ie){k^ - D) 
We now introduce the function g{A,C) defined as 



g{A,C) = ^ I dq 



Piq') 



nJo ' {q^-A){q^-C) • 
In terms of this function Eq. (8) can be written as 
Dm{A) = 1 - iA-D)NjesiD)giA,D) 
A-D rV^-§^fi^,(A,.^) 



(8) 



(9) 



(10) 



This a linear integral equation for Dj(s{A). Its linearity is 
an important property because it allows one to take more 
subtractions and still being amenable for an iterative solution. 
We take as a convenient subtraction point D = 0. In the case 



of the S waves, this choice relates the subtraction constant, 
Nj(s{0), to the corresponding scattering length, a,, by 



NjcsiO) = - 



4nas 



(11) 



The only elastic S wave is the '^o, and the scattering length for 
this wave is Us — —23 .758 ±0.04 fm. For this partial wave, we 
can write down the following integral equation: 



Dj(s{A) = \-ANm{Q)g{A,Q) + 



nj- 






(12) 



This is an integral equation for Dj(s{A) with A on the LHC, 
and it can be solved with the method exposed in the Appendix. 
Once Dj(s{A) is solved from Eq. (12), NjtsiA) is determined 
by inserting the former into Eq. (7), which now reads {D — 0) 



NMsiA)=Njes{0)- 



dk 

^ J-oo 



Ms{k^)Dm{k^ 
k^k^-A) 



(13) 



For a P wave, the same equations hold with A'yfs(O) = 
because in this case the amplitude vanishes at threshold as 
|p|2,andDj«(0) = l. 



C. Higher waves {( > 2) 

Equations (12) and (13), because of the full implementation 
of rescattering in Eq. (10), do not guarantee that the resulting 
partial wave amplitude behaves as A^ for A — > 0, with £^2. 
At LO Aji:s = AIL gives rise to OPE through the dispersive 



integral 



TJI^iA) = TjYsiO) + ^ £fM% ^ (14) 



ijl^iO) a subtraction constant. As discussed above. 



with Tj-i 

since the OPE amplitude [1] tends to constant for A ^• 
oo, the Sugawara and Kanazawa theorem requires that one 
subtraction is needed. The fact that for i > a partial wave 
vanishes as A^ for A ^ makes that TjY^{Q) = when £>0. 
This threshold behavior also implies that Ajf^ must fulfill the 
set of £ — 1 sum rules (constraints) 



dk 



,,2X 



0. 



(15) 



with A =2,3,...,^ and £ ^ 2. These constraints are obtained 
straightforwardly by performing the expansion of Eq. (14) in 
powers of A and imposing that Tjes{A) -^ A^ when A -^ 0. 

Let us now consider again Eq. (7). As Djfs{A) -^ 1 for 
A ^ then 7}fj(A) — > Nj(s{A) in this limit. The expression 
for A'jf5(A), Eq. (7), is similar to Eq. (14). Indeed, they would 
be the same equation if Dj(s {A) were replaced by 1 in Eq. (7) 
(and with Aj(s{k^) evaluated at LO). As a result, Njfs{A), 
determined by implementing Eq. (10) into Eq. (7), does not 
vanish as A^ for A — > 0, because of the departure from 1 of 
Dj(siA) in an actual calculation. 



It is convenient to proceed in a such a way that the 
right behavior of Tj(s{A) around threshold is incorporated 
explicitly. For that purpose we consider the N/D equation for 
Tjes{A)/A'^, instead of Tj(s{A). The quotient 7>«(A)/A* has 
no pole at A = because Tj(s{A)^A^ when A ^ 0. Notice 
also that due to unitarity Tj£s{A)/A^ — )- for A -^ +oo. Then, 
according to the Sugawara and Kanazawa theorem [32, 40], no 
subtractions are needed for the DR of Njts{A) with i>Q. EFT 
results do not always share the right high energy behavior so 
that subtractions will be certainly needed for a higher order 
calculation of Aj(s{A). At LO this is not the case because 
cxD. We then have (i > 0) 



Ajfs^ 1 /A for A 



7>«(A)-A 



e NjisjA) 
DjtsiA) ' 



1 r^ . AMs{k^)Dm{k^) 



(16) 
(17) 



A r<=o 



Dm{A)^\-- Aq 
It Jo 



= 1 + — / Ak- 

7t- J-oo 



,piq^)q^^'-'^Njts{q') 



Am{k^)Djes{k^) 



m 



X / dq' 



p{q)q 



2^^,2(t-\) 



[q^-A){q^-e-) 



(18) 



where the subtraction has been taken at threshold. 

The previous equation for Dj('s{A) is not satisfactory when 
£^2 because the last integration on the right-hand side (rh.s.) 
of Eq. (18) is divergent. In this way, by applying the N/D 
method to Tji!s{A)/A^ we have changed the problem of the 
bad behavior of Tj(s{A) around threshold into a high energy 
problem in the form of divergent integrals. To end up with 
a convergent DR for Djisi^) ™ ^l- (18) it is necessary that 
Njes{A) vanishes at least as 



Nj(siA) -^ l/A^' for A ^ OO 



(19) 



However, Njes{A) vanishes only as 1/A, independently of £, 
as follows from Eq. (17). The set of constraints needed to 
satisfy the asymptotic behavior in Eq. (19) can be deduced 
by performing in Eq. (17) a high energy expansion of Njts (A) 
in powers of 1/A. It results in 



L 

-CXD 



Ajisik^)Djts{k^ 

^2X 



0, 



(20) 



with X =2,3,... J, and £ ^ 2. These sum rules generalize the 
ones fulfilled by A]'^^{A) in Eq. (15).^ 

The usefulness of the £—\ restrictions in Eq. (20) can be 
well appreciated by rewriting Nj^si^) in Eq- (17) as 



1 



I: -2 



1 



Nm{A) = --Y. ^,„^, I 

m=0 



dt 



Ajes{k^)Dj(s{k^) 

l^2(e-m) 



nA^-^ 



L 



dk- 



, Ajes{k^)Dj(s{k^) 
k^{k^-A) 



(21) 



The last term on the rh.s. of the previous equation vanishes 
explicitly as 1/A^ for A — )• oo, while the terms in the sum on 
m are zero once the constraints of Eq. (20) are fulfilled. In this 
way, inserting this expression for A'7£5(A) in Eq. (18) one has 



Dm{A) = 1 + 



n J- 



V^z^£(^!l^-(^,(A,.^), 



k^ 



(22) 



and a convergent DR integral equation for Djiis{A) results. 

It should be stressed that Eqs. (16), (17), and (22) lead to 
the same equations as for the case of a f wave amplitude, 
£ = 1, cf. Eqs. (3), (12) and (13).^ In the case of a P wave, 
no constraints are needed because the right behavior near 
threshold is obtained straightforwardly. On the other hand, 
Eq. (22) can be readily applied to S wave by just adding 
the term proportional to Njes{0) present in Eq. (12). One 
subtraction should be taken in Eq. (17) in order to transform 
itasEq. (13)for^ = 0. 

Now, let us address the way to solve the N/D method, 
Eqs. (16), (17), and (22), in the presence of the constraints 
given in Eq. (20). It is well known [32, 37] that the function 
Dji:s{A) is determined modulo the addition of CDD poles 
[35]. These are associated to specific dynamical features of the 
interaction that arise independently of the LHC discontinuity, 
Ajcs{A), and unitarity. Typically, the addition of CDD poles 
corresponds to the existence of pre-existing resonances or 
to Adler zeros [31, 41]. Both facts are indeed absent in the 
low-energy A^A^ scattering [42]. We exploit this ambiguity in 
the function Dj(s{A) and include i— I CDD poles at infinity, 
so as to satisfy Eq. (20): 



Dj(s = 1 



n J- 



oo 1^ 



e-i 



f^^BiA-Bi 



(23) 



The last term in the rh.s. corresponds to adding the £ — \ 
CDD poles. The factor A/B,- in front of every CDD pole arises 
because the function Dj(s{A) is normalized to 1 for A = and 
it has the residue 7; at A = B,. The sum of the CDD poles gives 
rise to a rational fraction Q(_ 1 /f^_ 1 , where the subscript in Q 
and P indicate the degree of the polynomial in A. Since the 
only relevant fact at low energies is the ratio Ji/Bf we take at 
the end the limit B, -^ 00, with Ji/B^ not vanishing. The £—1 
CDD poles are gathered at the same point B and we write 



i-\ 



y, ^ AY!-J,cnA- 



BjA-Bi 



[A-B) 



(24) 



The coefficients c, are finally determined by requiring that the 
set of £ — 1 constraints in Eq. (20) are satisfied. The calculation 



^ Equation (19) is a consequence of Eq. (16) because for A — > +00, due 
to unitarity, the ratio Tjts{A)lA' tends to I/A'+'Z^ while Djts{A) -^ A'/^ 
(when only one subtraction is taken.) 



' It is equivalent to have the explicit factor A in Tj/s, Eq. (16), or included in 
the definition oiNj(s, Eq. (13). 



is performed in terms of finite but large B, and one has to 
check that the resuhs are stable by taking B arbitrarily large. 
At the level of low-energy A^A^ scattering we have modified 
Dji's{A) by adding a polynomial of degree i—l with fixed 
coefficients. 

We end with the following expressions: 



TjisiA) : 



' DMA) 



NjtsiA) = - 

K 



^ - Ans{k^)Djts{k^) 



DjisiA) = 1 

^Ln=oCnA 

{A -By-' 

and the constraints 



dk 

— OO 

A 

n 



i^jcsy 



k^^{k^~A) 
Ajfs{k^)D, 
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dk 

-OO 



i^jcsy 



JCS 



ik') 



k^ 



g{A.k^) 



(25) 



(26) 



(27) 



dk 



, AMs{k^)DMsik^) 

^2A 



0, A=2,3,...,£, e^2. (28) 



The previous formaUsm is also meaningful for the case in 
which Aji/s{A) — > C, with C a constant, as A ^ cxd. We 
do not discuss in this work its extension to the case when 
Aj(Js{A) diverges for A -^ oo, as we are interested now in 
LO A^A^ scattering. This generalization of our formalism will 
be discussed when considering higher orders in the chiral 
expansion of Aj(s{A) [43], which include the important TPE 
contributions [44, 45]. 

Summarizing the results of this section, we have presented 
a general approach based on the N/D method to construct 
A^A^ scattering partial wave amplitudes. For ^ = 0, 1 Eqs. (12) 
and (13) are employed, with Tj(js given by Eq. (3). For £ ^ 2, 
one has Eqs. (25)-(27), that must be solved in the presence 
of the constraints given in Eq. (28). The solution of this 
integral equation subjected to the constraints is discussed in 
the Appendix. 



III. RESULTS 

In this section we present our results for the phase shifts, 5, 
of the elastic partial waves with £ ^ 5. We compare them with 
the Nijmegen partial wave analysis (PWA) [42] in Figs. 2-4. 
Our results are represented by the solid (black) lines, and the 
Nijmegen data by the dash-dotted (red) lines. We show the 
results up to |p| = 300 MeV. Notice that the pion production 
threshold opens at |p| ~ 360 MeV and the three-momentum is 
no longer small, |p| ~ y/mm^ ^ m-,i. 

In Fig. 2 we show the lowest elastic waves, namely, '^o, 
^Pi, ^P\ , and ^Po> whose amplitudes do not contain CDD poles 
because (, <2. The agreement in the ^Pi and ^Pi partial waves 
is quite satisfactory. For the '^o it is known that a higher order 
chiral counterterm is needed to reproduce the large effective 
range and thus improve the agreement with the data [19]. In 
the case of the ^Pq wave, large corrections stem from TPE. 
Since this is a LO calculation, none of these corrections is 
included. However, an important point should be stressed: 
in this paper a regulator (cut-off) independent and unitary 
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FIG. 2. Comparison of the 'Sq, ^Po, ^P\, and ^Pi phase shifts with 
the Nijmegen PWA data. The solid (black) lines represent the results 
of this work, while the dash-dotted (red) ones represent the Nijmegen 
PWA [42], For the 'Sq and ^Pq, the dashed (black) lines correspond 
to the relativistic version of this work, see the text for details. 



description of the A^A^ interaction with the right analytical 
properties is reached. The agreement with the data can be 
improved by including higher orders in the LHC. 

For the '5o and ^Pq waves we have also tried with a 
relativistic calculation of the function g{A^k'^), Eq. (9), since 
these are the waves for which the discrepancies with the data 
are larger. In this approach, the p{q') function is replaced by 
its relativistic counterpart in the S'-matrix and in the integrals 
where it is involved. 



P{q^) = 



q^m 



4n 



Piq' 



'q^m 



A% 



V? 



The results obtained are represented by the dashed (black) 
lines in Fig. 2. Though the corrections are in the right 
direction, the discrepancies with the data are still large. As 
expected, relativistic corrections are small in the energy range 
shown, though noticeable for the ^Pq partial wave for |p| > 
200 MeV. 

For higher waves one needs to include £ — 1 CDD poles 
in order to fulfill the constraints in Eq. (28). Thus, it is 
guaranteed that the partial waves have the right behavior 
at threshold, vanishing as Ip^*^. Our results are shown in 
Figs. 3 and 4, and a good agreement with the data is 
achieved, except for the 'D2 wave. Our curves are quite 
similar compared with the LO results of Ref. [19]. This 
reference offers an approach with cut-off independent results 
with the NN potential (Vnn) given by OPE. The largest 
discrepancy concerns to the ^Pq partial wave where in 
Ref. [19] a counterterm is promoted from higher orders 
so as to achieve cut-off stable results for A ^- cxd due 
to the attractive l/r' tensor force in OPE. As a result. 
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FIG. 3. Comparison of the '£)2, ^D2, 'F3 and •'F3 phase shifts. For 
the notation, see Fig. 2. 
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FIG. 4. Comparison of the 'G4, ^Gix, ^H=, and ■'//g phase shifts. For 
the notation, see Fig. 2. 



the agreement with data is much improved. The main 
difference between our approach and that of Ref. [19] is 
the treatment of the LHC. Namely, for the ^Pq wave the 
iteration of the A^A^ potential is responsible for the need of 
this extra counterterm. The first iteration VnnGVnn (with G 
the unitary two-nucleon reducible loop function and V^n the 
nucleon potential) is a new source of LHC discontinuity [46] 
containing contributions from TPE and iterated OPE. The 
real part stemming from the former is divergent. Within our 
approach the sources of LHC discontinuity from VnnGVnn 
are NLO according with the standard chiral counting. At that 



order new subtractions are required [43] which will mimic the 
role of the extra counterterm taken in Ref. [19]. While our 
method is based on the calculation of Ayf5(A) perturbatively 
along the LHC, the application of a Lippmann-Schwinger 
equation with the chiral A^A^ potentials is based on the 
perturbative calculation of the latter [6-8]. In both cases 
the diagrams required for the calculation of Ay^^ and Va^a^ 
are two-nucleon irreducible, which justifies its perturbative 
treatment [1,4, 6-8]. In both cases as well the RHC is exactly 
resummed, as required because of the enhanced two-nucleon 
reducible diagrams. This resummation is performed in terms 
of the interaction kernel, Aj^s or Vnn, depending on the 
approach. The N /D method respects the LHC discontinuity 
so that Ajis is the same as in the final partial wave 
amplitude. For a Lippmann-Schwinger equation this is not 
the case as new sources of imaginary parts along the LHC 
result from the iteration of Vnn [46]. It is also worth 
stressing that our approach based on the N /D method is a 
dispersive one offering results that by construction are cut-off 
independent, while this is still an issue in the application of 
the Lippmann-Schwinger (or Schrodinger) equation to A^A^ 
scattering with Vm calculated from ChPT [19, 25, 27, 28]. 

For the '^o and 'D2 partial waves, for which we do not have 
good agreement at LO with data [42], our results are indeed 
very similar to those of Ref. [19], too. In the case of the ^P\ 
partial wave our phase shifts run closer to data at low energies 
than those of Ref. [19]. 

In order to show the independence of our results with 
the value of B, the position of the CDD poles, once this 
value is large enough, we show in Fig. 4 for the ^Ga, partial 
wave different lines corresponding to B = lQ"m\ for n = 
2,3,4,5, and 6. A narrow band is obtained despite the large 
variation in the values of B considered. 



IV. CONCLUSIONS 

We have applied the N /D method to A^A^ scattering from 
chiral perturbation theory. In this method the two cuts present 
in a A^A^ partial wave, the right-hand cut and left-hand 
cut, are separated in two functions, DjfsiA) and Njcs{A), 
with A the center-of-mass three-momentum squared. While 
Djf:s{A) has only right-hand cut the function Nj(s{A) has 
only left-hand cut. The A^A^ partial waves, Tjcs — Njcs/Dji/s 
(£ = 0, 1) and A^'Njts/Djis (^ ^ 2), are determined in terms 
of their discontinuity along the left-hand cut due to multipion 
exchanges, Aj^siA). At leading order, considered in this work, 
only OPE contributes. For D- and higher partial waves, 
with orbital angular momentum ^ ^ 2, one has to impose 
the proper behavior of a partial wave near threshold, such 
that it vanishes as A^ for A ^- 0. This gives rise to £ — 1 
constraints in the form of sum rules involving the functions 
Aj(s and Dji!^. Since the function Dj(s{A) is determined 
modulo the addition of Castillejo-Dalitz-Dyson poles (that 
correspond to zeros of the A^A^ partial waves along the real 
axis) we have added £ — 1 of such poles at infinity in Dji:^ 
for £ > 2. By sending such poles to infinity no zero at finite 
energies is included for any A^A^ partial wave. In addition, 
the residues of these poles in Dj^s ^rc fixed once the sum 



rules are satisfied, so that no new parameters are included. At 
low energies the Castillejo-Dalitz-Dyson poles behaves like 
adding a polynomial of degree £ — 2 to Z)yo- 

The resulting A^A^ partial waves do not contain any 
regulator A subtraction constant is required for the ^Sq 
partial-wave that is fixed by reproducing the experimental 
scattering length. Our results are very close to those of 
Nogga, Timmermans, and van Kolck [19] that provide cut-off 
independent A^A^ partial waves with OPE as potential. The 
only noticeable difference concerns the ^Pq partial wave for 
which Ref. [19] achieves cut-off stable results by promoting a 
higher-order counterterm to leading order due to the attractive 
1 /r^ tensor force in OPE. In our approach there is no special 
treatment for the ^Pq partial wave compared to others because 
of the perturbative treatment of the discontinuity across the 
LHC. Our results are a prediction for the ^Pq phase shifts at 
leading order For the ^Pi partial wave our phase shifts run 
closer to data than those of [19]. 

This method is ready to be extended to higher orders 
and it would be of great interest to apply it including 
TPE contributions to Aj^g. TPE is needed for an accurate 
description of the A^A^ data [44, 45]. Its extension to coupled 
channels with £ = Jztlis being studied [36]. 
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APPENDIX: SOLVING THE INTEGRAL EQUATION 

FOR D{k^) 

In this appendix, we focus on the solution of the integral 
equation (27) subject to the constraints (28). For simplifying 
the discussion we drop the subscripts J£S. The integral 
equation Eq. (27) can be written in a compact way as 

D(A) = 1 + - f'dk^^^l^)^g(A,k^)+hiA), (Al) 

where g{A,k^) is defined in Eq. (9) and 



(-2 



h{A) = 



(a^bFT,?/'^' 



(A2) 



Let us introduce the function d{A) as D{A) = d{A) +h{A), 
that is, d{A) is the piece of D{A) that does not contain the 
CDD poles. As a first step, we write the coefficients c, in 



terms of the d{A) function. This is done by writing the sum 
rule constraints Eq. (28) in terms of d{A), giving rise to 



/,= 



'^^^^, A{k^)d{k^) ^ f^^^^,A{k^)h{k^ 



— oo 

-2 



k^k^' 



k^k^> 



dk' 



A{k^)k^j 



k^k^'{k^--B] 



(-1 



(A3) 



(A4) 



where we have shifted A = / + 2, so that / runs from to 
£ — 2. Note that the integrals niij can be calculated directly 
for a given A(A) in terms of B because the unknown function 
d{A) does not appear in their calculation. Thus, 



1-2 



7=0 






(A5) 



being m the inverse of the matrix ||ot,j||. Next, we rewrite 
the integral equation Eq. (Al) in terms of d{k^), and insert 
Eq. (A5) for the coefficients c,. It results 



ij=0 



dq 



^ d{q^)A{q^ 



(A6) 



Now, by interchanging the dummy integration variables k^ 
and q^, we can finally write 



<,W = i + ^/ 






-8{A,k^)) 



(A7) 



f-2 



8{A,k^ 



.k^ 



mJ l%2H,')^(^,,') 



U.0- -- {q'-BY-^ 



We have now an integral equation for the d{k^) function that 
depends on known functions. It can be written in a compact 
way as 



^(A) = l+ / dk^f{A,k^)d{k^) 



(A8) 



It is convenient to perform a change of integration variable 
so that one ends with a finite integration domain, e.g. with 
X — l/k^. In this way 



d{A) = 1 + / dxf{A,x)d{k^{x)) 



(A9) 



■Tl 



This is an inhomogeneous Fredholm integral equation. We 
solve it numerically by discretizing the integral on it, 

d{Ai) - 1 +£/(A;,X,)«(^7M(fc^(^7)) . (AlO) 



where the (o{x) function is the weighting function taken for 
the integration. By calUng d{k^{xi)) = dt, f{Ai,Xj)(o{xj) = 
rjij, this equation can be recast as 



L(^'7 - ^-Mj =- 1 , 



(All) 



which is a linear equation, that can be solved by standard 
methods, giving the desired function d{A). To obtain D{A), 
the function h{A) must be added, but this is also a direct task, 
since the c,- coefficients can be calculated once d{A) is known, 
Eq. (A5). Of course, if no constraints must be satisfied, as it 
is the case for S and P waves, the same formalism with the 
corresponding simplifications should be used. 
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